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ABSTRACT 

This  paper  deals  with  the  numerical  solution  of  eigenvalues  of  a  three -parameter  Sturm-Liouville  problem  using 
Moment  method.  The  resulting  algebraic  equations  obtained  in  applying  Moment  method  on  the  problem  are  solved  to  find 
the  rough  estimates  of  the  eigenvalues  of  the  problem.  Rough  estimates  are  used  as  starting  approximations  in  the 
corresponding  shooting  method  to  obtain  their  actual  values. 
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1.  INTRODUCTION 

Multiparameter  eigenvalue  problems  (MPEVP's)  are  generalisation  of  one-parameter  eigenvalue  problems. 
Separation  of  variables  for  a  partial  differential  equation  leads  to  a  multiparameter  systems  of  ordinary  differential 
equatioins  (Roach(1974)). 

It  has  been  observed  that  one-parameter  problem  have  gained  much  development  both  theoretically  and 
numerically  in  comparison  to  multiparameter  problem.  Unlike  theoretical  aspects  (Atkinson  (1968),  Collatz  (1968), 
Fox  et  al.  (1972),  Roach  (1974),  Sleeman  (1972)),  the  numerical  methods  (Baruah(1987),  Baruah  (2011))  developed  for 
solving  multiparameter  problems  are  very  limited.  In  this  paper,  we  will  investigate  a  three-parameter  eigenvalue  problem 
in  the  form  of  associated  with  four  point  boundary  conditions. 

ybd  +  U+  ff/Cs)  +  vgidhte)  =  0  (1) 
y(a>  =  y(b>  =  y(c)  =  yid)  =  0  (2) 

Where  a  <  b  <  C  <  d,  f(x}  and  &(x)  are  given  real  valued  continuous  functions  of  the  independent  variable 
x  E  [q,  ri].  The  values  of  the  parameters  A,  u  and  v  are  called  eigenvalues  and  the  non-trivial  solutions  y(x}  are  called 
eigenfunctions  of  (l)-(2)  corresponding  to  the  eigenvalues. 

The  parameters  are  to  be  so  chosen  that  the  above  equation  should  have  non-trivial  solution  satisfying  boundary 
conditions,  and  also  the  parameters  have  prescribed  number  of  zeros  in  each  of  the  intervals  {b,c)  and  (c,  d} 

(Sleeman  (1972)).  The  points  at  which  y(x}  vanishes  in  [n  b],  [h,  c]  and  \c,  d]  are  called  zeros  of  the  eigenfunction.  In  this 
paper  rough  estimates  of  the  eigenvalues  are  obtained  by  an  application  of  Moment  method  which  is  based  on  the  principle 
of  Weighted-Residual  methods.  Using  the  rough  estimates  as  starting  approximations  in  the  corresponding  shooting 
method  of  the  problem,  actual  values  are  obtained. 
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2.  FORMULATION  OF  THE  METHOD 

Initially  the  complex  partial  differential  equations  were  solved  by  finite  difference  method.  But  it  is  difficult  to 
use  finite  difference  method  for  irregular  geometries.  To  overcome  this  difficulty  a  better  method  was  developed  which  is 
known  as  finite  element  element  method.  The  concept  of  FEM  is  based  on  the  classical  variational  method.  But  there  are 
some  limitations  on  classical  theory  to  apply  on  more  types  of  problem  while  the  Weighted-Residual  methods  provides  a 
strong  basis  for  wider  range  of  problems. 

In  FEM  value  of  some  kind  of  integral  is  minimised  over  the  entire  domain  and  the  process  of  minimisation  is 
carried  out  element  wise  to  get  an  over  all  minimum  (Baruah(201 1)).  Moment  method  which  belongs  to  Weighted 
Residual  methods  is  also  one  such  approach  in  which  we  seek  some  trial  function  to  obtain  an  approximate  solution. 
Therefore,  we  will  give  a  brief  account  Weighted  Residual  method  followed  by  Moment  method  of  our  present  discussion. 

3.  WEIGHTED  RESIDUAL  METHOD 

The  principal  idea  of  weighted  residual  methods  is  to  find  an  approximate  solution  to  differential  equation  in  the 
form  (M.  K.  Jain  (1983)  and  K.  S.  Rao  (2007)), 

JrW  =  /  (3) 

Subject  to  the  boundary  conditions  such  as 

£/u  =  g}  (4) 
In  order  to  find  an  approximate  solution  to  the  above  boundary  value  problem,  we  introduce  a  set  of  functions  of 

the  form 

w  =  M[j  +  £jUGjfi  (5) 

where,  ijfij-  are  the  linearly  independent  functions.  Also  4tp  satisfy  inhomogeneous  boundary  conditions  and  cjbj- 
satisfy  homogeneous  boundary  conditions.  The  coefficients  Cj  are  unknowns  which  are  to  be  determined  by  solving  a 
system  of  equations.  Since  (5)  is  an  approximate  function,  we  observed  that,  when  substituted  in  (3),  it  will  not  satisfy 
exactly.  Thus,  we  get  an  error  or  equation  residual  denoted  by  F(jf,  b]  and  can  be  written  as, 

ti&-f  =  Efra)  (6) 

This  gives  the  measure  of  the  extent  to  which  the  function  F'ix,  d)  satisfies  the  differential  equation.  As  the 
number  n  of  the  functions  <p  .-(x)  is  increased  in  successive  approximations,  we  hope  that  the  residue  E(x,  Cl]  will  become 
smaller.  The  exact  solution  is  obtained  when  the  residue  is  identically  zero.  Since,  it  is  not  possible  to  make  the  error 
E(x,  a)  identically  zero,  we  shall  try  to  make  E(x,  a}  as  small  as  possible  in  some  sense.  In  the  error  distributing  methods, 
the  error  function  E'ix..  a)  is  made  to  approximate  the  zero  function,  in  the  sense  that  each  of  the  weighted  integral  of  the 
residue  of  E(x,  c0  with  respect  to  the  weight  functions  Wp  j  =  L  2, ... ,  ft  is  equal  to  zero. 

jE(x,  a)Wjdx  (7) 

The  approximate  methods  based  on  (7)  are  called  the  weighted  residual  methods.  There  are  many  ways  of 
choosing  the  weighting  functions       Their  choices  lead  to  different  methods.  They  are 
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•  Galerkin  method 

•  Least  square  method 

•  Method  of  moments 

•  Collocation  method 

As  already  stated,  in  this  paper  we  consider  Moment  method  based  on  this  principle  as  described  below. 
3.1  Moment  Method 

In    this    method,    the    weighting    functions    are    chosen    from    family    of    polynomials    such  as 
wf=i;,    i  =  1,  2, .,,  ..n  —  1  so  that 

_f  Eix,  a}x^x  =  0 ,   £  =  1, 2, ,  „ ,  n  - 1  (8) 

This  gives  a  system  of  linear  or  non-linear  equations  for  the  solution  of  the  parameters  D]_,  q:,  ,,, ,  oR. 

4.  NUMERICAL  EXAMPLE 

The  numerical  example  of  the  problem  (l)-(2)  considered  here  is  a  three-parameter  problem  defined  by  the 
equation  (Collatz  (1968)), 

y"bd  +  (1  +  2ff  -cos  x+  Zv  cos  Zx)y(x)  =  0  (9) 
Subject  to  the  boundary  conditions 

ytO)  =  yW  =  y{2)  =  y(3}  =  0  (10) 
Here  we  choose  the  approximate  solution  of  (9)  satisfying  the  boundary  condition  (10)  as 

kGt)  =  a.lX(3  -  x)tx  -  1>(2  -x)  +  o:x:(3  -x)Cc  -1>(2  -  x)  4  o,**(3  -*•)(*- l)  (2  -  x)  (11) 
Substituting  the  approximate  solution  (1 1)  in  the  problem  (9)-(10),  we  get, 

Htjf.a)  =  a^ClZx1  -  '3bx  +  W  +  Tlx'  +  bbx  -  12)  +  at(iQx*-T/.9x*  4  132xa  -  Mix)  +  (A  + 

2a  cdsx  +  2vccs  Zx){%lx4  -  6x*  +  Ux*  -  fist)  +  GjU"  -  6x*  +  lis*  -  fix2)  +         -  fix*  +  11**  - 

6*^}  (12) 

For  exemplification  the  weight  functions  are  considered  as  \V^_  =  1,  IVj  =  X  and  IVj  =  x2. 

Using  equation  (8)  in  the  intervals  [0.,lL  [1^2]  and  [2.,  3]  we  obtain 

SgEtx,al  WLdx  -  (-0.63333A  -1.1150*,u  -  0.72177v  +  S.000O0)aL  4  (-0.2S333A  -  0 .47185*  - 
0,231&2v  4  2,0000  0)  a-  4  (-0,157m  -Q,249S3,u  -  0,1351  6v  +  2,00000)  b3  =0  (13) 

/^(x  a)  H^iv  -  (0,266661  4  0,05062  a  -  0,6S772v -4,D0000)Ol  +  (0  .E500Q1  4  0 ,040S3ii  -  0,9959  Qv  - 
&.JUUUUJD;  +  (L.^tM  4-  U-Ull^aa  -  l,53Kttv  -  10,UULU0)sg  =0  (14) 

l:!F(xs]  WLdx  =  (-0,633331  4 1,02845  u  -  0,4431  Ov  4  3,00000}  oL  4  (-1,616671-  4  2.65950,ii  - 
l,23744v  +  22,00000)02  +  (-4,157141  4  6.92259  u  -  3.44479-v  4  62.00000)q,  =  0  (15) 

^E(x,a}Wtdx  =  (-0.23333,1  -  0,47135 u  -  0,23162v  4  2,00000}  oL  4  (-0,157141  -  0,24933 a  - 
0,03969v  4  2,00000)  B;  4  (-0,098311  -0.15116u  -  0.03757V  4  2,00000)d,  =  0  (16) 
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u>IVaiii  -(0. 5 SO 00 A +  0.040 9 3u  -  0,93593l'  -&.0000D)bl  +  (0.342S01  +  k01129p  -  l,33SZ8if  - 
IQ.OQOCOJOj  +  U.317St>l-  D.05S32U  -  Z.4UJ9V  - 13,0000 0}q2  =0  (17) 

J^E(x,  a}Wzdx  =  (-1. 61667 Jl  +  2.65950.U  -  1.23744V  +  22.0000O)a1.  +  (-4.157141  +  6.9225%  - 
3,44479V +  62.00000)as  + (-10.765431  +  13.13230  u  -  9.57093v  +  173.00000] o3  =0  (18) 


^F(x,  a}  W3dx  =  (-0.157141  -  0.24933 u  -  0,08969v  +  0.73333)  gl  +  (-0.093311  -  0.15116.U  - 
0.037571'  +  1.4333 3] a,  +  (-0.067461  -0.09993.U  -  0.0155  6v  +  1.63571) o3  =  0 


(19) 


^E(^a-}W,dx  =  CQ.Qi2B6A  +  0.01129a  -  1.5362Bv  - 9 ,26S67)cL  +  (1.317S61  -  0.058-32.41  -  2.41119v  - 

l&,90O0O)t3n  +  (2.09921A  -  0,205 73«  -  3.8461+v  -32,31429)q3  =0  (20) 

S?E(x,  a}  Wtdx  =  (-4.157141  +  6.92259a  -  3.444S0v  -  60.733333)  aL  +  (-10  765431  +  18.13230j*  - 
9,57093v+  174,76667}     +  (-23,067*61  +  47, 7S1 3  4ft  -  26.56*3 lv  +  509, 63570)^  -  0  (21) 

Eliminating  the  parameters  BL,  o:  and  c.s  from  the  equations  (13)-(21)  taking  three  at  a  time  we  get  different 
algebraic  equations  of  fourth  order.The  algebraic  equations  are  solved  (S.  S.  Sastry  (2006))  to  get  the  starting  values  of  the 
eigenvalues  (1ft  v)  of  the  problem  (9)-(10)  as  shown  in  the  table-1  below. 


Table  1:  Starting  Appox  for  Shooting  Method 


Equations  Considered  to 
Eliminate  a,,  a  -  and  a-. 

Annrovimntr  Solutions:  of* 

(A,  fi/  v  j 

(13),  (14),  (15) 
(16)  (17)  (18) 
(19),  (20),  (21) 

(0.39355,  (-0.17528,  7.66561) 

(13),  (20),  (21) 

(13)  ,  (15),  (17) 

(14)  ,  (15),  (16) 

(0.83488,  -0.32253,  3.57482) 

(14)  ,  (15),  (16) 

(15)  ,  (16),  (17) 

(16)  ,  (17),  (18) 

(0.06837,  0.05869,  2.43435) 

(13),  (14),  (15) 
(13),  (16),  (17) 
(13),  (14),  (16) 

(2.93326,  -0.08234,  19.90332) 

(13)  ,  (18),  (19) 
(15),  (17),  (21) 

(14)  ,  (16),  (20) 

(-1.58854,  -0.22714,  4.17804) 

(13)  ,  (16),  (21) 

(14)  ,  (17),  (19) 

(15)  ,  (18),  (20) 

(0.22315,0.02355,8.76451) 

(13),  (17),  (21) 
(15),  (17),  (19) 
(15),  (16),  (18) 

(1.65694,  -0.01337,  -2.05649) 

(16),  (20),  (21) 
(15),  (19),  (20) 
(14),  (18),  (19) 

(1.63134,  1.09124,  -4.53315) 

(14),  (15),  (19) 
(18),  (19),  (20) 
(13),  (16),  (21) 

(-0.68440,  -0.56030,  16.81688) 

(13)  ,  (17),  (21) 
(16),  (19),  (21) 

(14)  ,  (17),(21) 

(0.31307,  -0.00003,  8.58012) 

(13),  (17),  (18) 
(16),  (20),  (21) 
(15),  (17),  (19) 

(1.03074,-1.66158,4.55446) 

(13)  ,  (16),  (19) 

(14)  ,  (17),  (20) 

(15)  ,  (17),  (21) 

(0.10510,  5.64528,  -12.95556) 
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The  rough  estimates  from  Table  1  are  used  as  starting  approximation  in  the  corresponding  shooting  metrhod  to 
obtain  the  actual  values  of  the  eigenvalues  (A,f£,v"i  of  the  problem  (l)-(2).  The  distinct  actual  eigenvalues  are  found  as 
(9.85032,  1.17066,  0.21466),  (20.32325,  -7.60925,  4.53812),  (37.99081,  -22.74049,  12.21601), 
(49.93017,  -1.15909,  -29.31049),  (59.66886,  -14.81822,  -21.84090),  (84.06648,  -34.78913,  -11.11194),  the  corresponding 
eigenfunctions  have  number  of  zeros  in  the  intervals  (0,  1),  (1,  2)  and  (2,  3)  are  (0,  0,  0),  (0,  0,  1),  (0,  0,  3),  (1,  3,  1), 
(1,  3,  3),  (1,  3,  5)  respectively. 

In  the  numerical  example  illustrated  above  the  starting  approximations  of  the  eigenvalues  depend  on  the  choice  of 
approximate  solution.  The  problem  defined  in  the  domain  [a,  ti]  are  considered  separately  in  the  intervals  [a,  b],  [b,  c]  and 
\c,  d]  so  that  enough  algebraic  equations  can  be  constructed  to  obtain  the  unknowns  a  "s  together  with  A,  u  and  v. 

4.  CONCLUSIONS 

Compared  to  theoretical  aspects,  numerical  treatment  of  three-parameter  eigenvalue  problems  are  still  very  few 
and  hence  some  contribution  in  this  area  are  always  in  need.  In  this  paper  eigenvalues  of  a  three-parameter  Sturm-Liouville 
problem  is  computed  using  Moment  method,  which  fall  under  Weighted-Residual  method.  One  of  the  most  important 
advantage  of  this  method  is  that  it  can  be  applied  to  both  non-selfadjoint  and  non-linear  problems(Finalayson  (1966)). 
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